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The influence of antiferromagnetic order on the superconductivity in the non- 
centrosymmetric heavy fermion compound CePtaSi and related materials is discussed. Based 
on our RPA analysis for the extended Hubbard model two phases could be stabilized by a 
spin fluctuation induced pairing, with either dominantly p-wave or d-wave symmetry. The an- 
tiferromagnetic order plays an essential role for the low-energy physics, in particular, for the 
appearance of line nodes in the gap and the enhancement of spin susceptibility below T c . Various 
properties and possible phase diagrams under pressure are analyzed. The present experimental 
situation suggests that the p-wave phase is most likely realized in CePtaSi. 

KEYWORDS: Superconductivity without inversion center; antiferromagnetic superconductor 
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Since the discovery of superconductivity in the non- 
centrosymmetric heavy Fermion compound CePtaSi, 1 
superconductivity in materials without inversion center 
has been attracting growing interest. Many new non- 
centrosymmetric superconductors with unusual proper- 
ties have been identified among heavy fermion systems 
such as Ulr, 2 CeRhSi 3 , 3 CeIrSi 3 , 4 CeCoGe 3 5 and others 
like Li2Pd a; Pt3_ a ;B, 6 and KOs20g. 7 One immediate con- 
sequence of non-centrosymmetricity is the necessity for a 
revised classification scheme of Cooper pairing states, as 
parity is not available as a distinguishing symmetry. The 
pairing states is considered as mixtures of states with dif- 
ferent parity imposed by the presence of antisymmetric 
spin-orbit coupling (ASOC). 8 Recent theoretical studies 
let to the proposal of various intriguing properties of such 
a superconductor. 8-14 

In the past the relation between superconductivity 
and magnetism has been one of the aspects of ma- 
jor interest in heavy fermion systems. Interestingly, all 
presently known non-centrosymmetric heavy Fermion su- 
perconductors, i.e. CePtaSi, Ulr, CeRhSi3, CeIrSi3 and 
CeCoGe3, coexist with the magnetism. Although mag- 
netism affects the electronic state profoundly, most of 
the theoretical studies except for Refs. 14 and 15 ne- 
glected this aspect so far. The aim of the present study 
is to elucidate how the magnetism influences the su- 
perconducting (SC) phase and how it may be involved 
in deciding the pairing symmetry in CePtaSi. Among 
the non-centrosymmetric heavy fermion superconduc- 
tors, CePtaSi has been investigated in most detail be- 
cause the superconductivity occurs at ambient pressure. 1 
Although we focus here on CePtaSi, we believe that some 
of our results are qualitatively valid for the other com- 
pounds too. 

In CePt 3 Si superconductivity with T c = 0.75K ap- 
pears in the antiferromagnetic (AFM) state with Neel 
temperature Tn = 2.2K. 1 Neutron scattering measure- 
ments characterize the AFM order with an ordering 
wave vector Q = (0, 0, n) and magnetic moments in 



the o6-plane of the tetragonal crystal lattice. 16 The na- 
ture of the SC phase has been characterized by sev- 
eral experiments. The low-temperature properties of the 
thermal conductivity, 17 superfluid density 18 and specific 
heat 19 indicate line nodes in the gap, while the coher- 
ence peak in NMR \jT\T is a feature expected rather 
for a conventional superconductor. 20 The upper criti- 
cal field H C 2 ~ 4T exceeds the standard paramagnetic 
limit, 1 which seems to be consistent with the Knight shift 
data displaying no decrease of the spin susceptibility be- 
low T c for any field direction. 21 ' 22 The combination of 
all these features is incompatible with the usual pairing 
states such as the s-wave, p-wave or cZ-wave state, and 
calls for an extension of the standard working scheme. 

For the following study of superconductivity in 
CePtaSi, we introduce the single-orbital Hubbard model 
including AFM order and ASOC, expressed as 
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U2_^n lA n ld . (1) 



We consider a simple tetragonal lattice and assume the 
dispersion relation as, 

e(k) = 2ti(cos fc x + cos fc y ) + 4^ cos fc x cos k y 

+2£a(cos 2fc x + cos2fc y ) + [2<4 + 4i5(cosfc x + cos fcy) 

+4^6(cos 2fc x + cos 2fc y )] cos k z + 2t*t cos 2fc z — \i, (2) 

which reproduces the so-called /3-band of CePtaSi 
as obtained from band structure calculations with- 
out taking AFM order into account. 23-25 The (3- 
band has a substantial Ce 4/-electron character 25 
and the largest density of states (DOS), namely 70% 
of the total DOS. 23 We determine the chemical po- 
tential fi so that the electron density per site is 
n and the parameters as (fi, £2, *3, *4, *5, ^ ^ = 
(1, -0.15, -0.5, -0.3, -0.1, -0.09, -0.2, 1.75) defining t x 
as the unit energy. 
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The second term in eq. (1) describes the ASOC due to 
the lack of inversion symmetry. In case of CePtaSi the g- 
vector has the Rashba type structure. 9, 10 Although the 
detailed momentum dependence of the g-vector is not 
easily obtained from band structure calculations, it can 
reasonably be expressed in terms of velocities v^ y {k) = 
ds(k)/dk x , y : g(k) = (—v y (k),v x (k),0)/v. We normalize 
g by the average velocity v [v 2 = J2k v x(fc) 2 + v y (k) 2 }. 
This form reproduces the correct symmetry and period- 
icity of the g(k) within the Brillouin zone. We choose the 
coupling constant a = 0.3 so that a band splitting due 
to ASOC is consistent with the band structure calcula- 
tions. 23 Figure 1 shows the Fermi surfaces in our model. 

The AFM order enters in our model through the stag- 
gered field Hq without discussing its microscopic origin. 
The phase diagram under pressure implies that the AFM 
order is mainly carried by localized Ce 4/-electrons which 
have a character different from the SC quasiparticles. 
The SC T c is little affected by the AFM order which 
vanishes at P ~ 0.6GPa 26 in contrast to the other Ce- 
based superconductor. 27 The experimentally determined 
order corresponds to Kq pointing along the ^-direction 
with a wave vector Q — (0, 0,7r). 16 For the magnitude 
we choose /iq -C W where W is the band width since 
the observed moment ~ 0.16^b is considerably less than 
the full moment of the Ce-ion. 16 We do not touch the 
complex heavy Fermion aspect, i.e. the hybridization of 
the conduction electrons with the Ce 4/-clcctrons form- 
ing the strongly rcnormalizcd quasiparticles. However we 
consider the Hubbard model as a valid effecive model to 
describe the low-energy quasiparticles. 28 

The undressed Greens function for U = has the ma- 
trix form, G(k,iu) n ) — (ko n l — H(k)) , where 



G(k,iuj n ) = 



H(k) 



G\k,iuj n ) G 2 (fc, iuj n ) 
G 2 {k +1 \Lo n ) G 1 (fc + ,iw„) 

e(k) -/i Q <7« \ 

-h Q aW e(k+) i ' 



(4) 



with e(k) = e(k)a ( °) +ag(k)a and k + = k+Q. G*(fc, iw n ) 
is a 2 x 2 matrix in spin space, w n = (2n + 1)ttT and T 
is the temperature. 

We turn to the SC instability which we assume to arise 
through electron-electron interaction incorporated in the 
effective on-site repulsion U. The linearized Eliashberg 
equation is obtained in the standard procedure: 

^^p,s 1 ,s 2 (k) = — ^ r p,q-s 1 ,s 2 ,s 3 ,s i (k, k')Tp q ^ S3 ^ S/L (k'\§) 

k',q,S3,S4 

' l l J P,Sl.S2 (fc) — y ' ( f > p,i.j,S 1 ,S3,S 2 .S i 

(k)A q (k"), (6) 

i,j,S3,Si 

where q = p (q = 3—p) for i = j (i ^ j), k" = k + (i — l)Q 
and 

n 

xGi 3 , S4 (-k+( P -l)Q,-iuj n ) (p=l,2). (7) 



Here, we adopt the so-called weak coupling theory of su- 
perconductivity and ignore self-energy corrections and 
the frequency dependence of effective interaction. 28, 29 
This simplification affects the resulting transition tem- 
perature but hardly the symmetry of pairing. 28 The effec- 
tive interaction V^ 9iSljS2 , S3jS4 (fc, k') originates from spin 
fluctuations which we describe within the RPA. 29 

The linearized Eliashberg equation allows us to deter- 
mine the form of the leading pairing instability which is 
attained for the temperature at which the largest eigen- 
value A in eq. (5) reaches unity. We perform the cal- 
culation at a given temperature, T = 0.02 and deter- 
mine the most stable pairing state as the eigenfunction 
of the largest eigenvalue for the sake of numerical accu- 
racy. 28 The typical value of the eigenvalue at T = 0.02 
and C/ = 4 lies around A = 0.4 — 0.6. 
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Fig. 1. The Fermi surfaces of the Hubbard model (eq. (1)) at 
a = 0.3 and /iq = 0. The cross sections at k z = ^, k z = ^ and 
fc z = are shown from the left to the right. 



First of all, we discuss the symmetry of the SC state 
which we express by the following extended parameteri- 
zation of the gap function: 



(3) A w (fc) = 



-d K (k) +id y (k) 



$(k) + d z (k) 

dx(k)+id y (k) 



(8) 



where we use the usual even parity scalar function $(fc) 
and the odd parity vector d(k). In the presence of AFM 
order a second component A 2 ,;s, S '(A:) appears. However, 
the basic properties and symmetries are little affected 
by A 2iSjS '(/c). Within the described scheme we identify 
two stable solutions of the Eliashberg equation eqs. (5-7). 
One pairing state has dominant p-wave symmetry whose 
order parameter has the leading odd parity component 
d(k) = (— sinfc y ,/3sinA: x ,0) and the admixed even parity 
part <&(fc) = cos k x + cos k y . The parameter [3 is unity in 
the absence of AFM order. 

The other stable solution has dominantly <i-wave 
character and can be viewed as an inter-layer pairing 
state: = {sin k x sin k z , sin k y sinfc z } (two- fold de- 

generate) admixed with odd-parity component d(k) ~ 
$(fc)(sin fc y , sin£; x , 0). In the paramagnetic phase the 
most stable combination of the two degenerate states is 
chiral: $±(fc) = (sin/c x ± isinfc y ) sinfc z which gains the 
maximal condensation energy in the weak-coupling ap- 
proach. In the AFM state, however, the two states of 
$(fc) are no longer degenerate. 

A brief view on the pairing mechanism clarifies the 
origin of the stable pairing states. The static spin sus- 
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ceptibility is peaked around Q = (0, 0, tt) consistent with 
the AFM order of the Ce moments. The in-plane fer- 
romagnetic correlations favor the in-plane p-wave pair- 
ing. On the other hand, the interlayer AFM correlation 
drives interlayer spin singlet pairing. The stability of the 
two states is determined by the band filling n and the 
Coulomb repulsion U . Small U and large n favor the p- 
wave state, while the d-wave state is more stable for large 
U and small n. The two states are essentially degenerate 
for U — 4 and n = 1.75. 




Fig. 2. (Color online) DOS p(e) for U = 4 in (a) dominantly p- 
wave state and (b) dominantly d-wave state. We show the results 
for e > because p(e) is particle-hole symmetric owing to its 
definition. 



Next we turn to the influence of the AFM order on 
these pairing states. As a first point we discuss the quasi- 
particle DOS which is obtained by diagonalizing the 8 x 
8 matrix, 



H s (k) 



H(k) 

-At(jfe) 



where 



A(fc) = 




-A(fc) 

A 2 , s , s '(fc) 
A MiS /(fc+) 



(9) 



(10) 



The DOS is obtained by the eigenvalues as p(e) — 
Tn Si k fi( £ ~ Ei(k)). We determine the momentum and 
spin dependences of order parameter within the lin- 
earized Eliashberg equation at T = 0.02 and assume 
that those structures do not change below T c . For our 
purpose it is not necessary to calculate the magnitude 
of the gap functions self-consistently as we are mainly 
interested on qualitative properties arising from the gap 
structure. Thus we choose the magnitude of the maximal 
gap, A g = 0.1, which may be large compared to the en- 
ergy scales a or Iiq. However, we adopt this value for the 



sake of numerical accuracy, having confirmed that lower 
values of A g do not alter the result in a qualitative way. 

We first consider the p-wave state (Fig. 2(a)). In the 
absence of ASOC and AFM order there are only point 
nodes along the [001]-direction leading to a quadratic 
energy dependence of the DOS: p(e) — c\S 2 . The in- 
clusion of ASOC (a = 0.3) yields two kinds of the 
line node. The admixture of the s-wave component is 
one cause of line nodes as discussed by Frigeri et al. 10 
The other origin of nodes lies in the specific structure 
of the g-vector. For the assumed band structure and g- 
vector, singularities of g appear not only along the [001]- 
direction (given by symmetry) but also accidentally on 
lines around (k x ,k y ) = (±0.47r, ±0.47r). The SC p-wave 
gap vanishes along the lines where d(k) _L g(k). 30 This 
second type of line nodes, however, depends strongly on 
details of material parameters. Within our model the 
length of the line nodes arising from these mechanism 
are short leading only to a weak linear energy depen- 
dence of the DOS, p{e) = ci£ as can be seen in Fig. 2(a). 
The linear DOS increases remarkably through the AFM 
order (Fig. 2(a)), caused by two effects: (I) the Bril- 
louin zone folding at k z = ±tt/2; (II) the modification 
of SC order parameter. Effect (I) has been investigated 
by Fujimoto. 14 It turns out that this effect is of minor 
quantitative importance, if /iq <C W as in the present 
situation. Hence the main effect originates from (II), be- 
cause the d x - and d y -components are no longer equiva- 
lent in the AFM state. The anisotropy parameter (3 in 
d(k) = (— sin/c y , /3sinfc x , 0) decreases with growing Iiq. 
In this way the SC gap becomes anisotropic leading to 
the extension of the line nodes and an even more pro- 
nounced linear energy dependence of the DOS. It should 
be noted that the line nodes discussed for this case are 
not symmetry protected but " accidental" . 

The d-wave state has already line nodes for symme- 
try reasons. In this case the low-energy DOS is not so 
strongly affected by the ASOC as shown in Fig. 2(b). 
But the slope increases, if the AFM order is included, 
since the pairing state changes its form from d xz ± id yz 
to d xz . The latter has obviously more line nodes. 




Fig. 3. (Color online) Schematic phase diagram in the P-T plane, 
(a) p-wave and (b) <i-wave state. "D" ("cD") shows the dxz _ wave 
(dxziidyz-wave) state. The critical temperatures of SC and AFM 
orders are written so as to be consistent with the experiment. 26 



Consequently, the low energy excitations are increased 
by the AFM order for both the p-wave and d-wave states. 
The resulting line node behavior is consistent with the 
experimental results at ambient pressure. 17-19 
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The p-wave case leads to a simple phase diagram in 
the P-T plane (Fig. 3(a)). However, the situation is more 
intriguing for the d-wave case, since there is an additional 
phase transition line meeting at the crossing point of 
T/v and T c . Generally we would expect an additional SC 
phase transition within the SC phase leading to chiral 
(i-wave phase at low enough temperatures, while the SC 
high-temperature phase has only a finite d xz -component 
(Fig. 3(b)). Although a second SC transition has been 
observed, 31 it remains unclear whether it represents an 
intrinsic property or is caused by sample inhomogeneity. 
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Fig. 4. (Color online) Magnetic susceptibility along the a- 
(dashed) and fe-axis (solid) against the staggard spin polariza- 
tion in = | < ^2 S s , a j ct s c is i > |. The thick and thin line 
show the results in the normal state and SC state at T = 0, 
respectively. 



eluding the ASOC and AFM order. For this purpose 
we chose the /3-band of CePtsSi, and found based on 
spin fluctuation mediated pairing interaction that the 
in-plane p-wave and inter-plane ci-wave states are most 
likely candidates. Both states show line node behavior, 
consistent with experiments at ambient pressure. 17-19 
The AFM order plays an important role in various re- 
spects: (I) the SC gap structure is remarkably deformed 
through the AFM staggered moment. (II) The AFM or- 
der can give rise to multiple SC phase transitions. (Ill) 
The in-plane magnetic susceptibility in the SC state can 
be increased giving an explanation for the Knight shift 
measurements. 21,22 The p-wave state is more likely re- 
alized in CePtaSi because it can explain the coherence 
peak in the NMR 1/TiT. 11 - 14,20 Finally, our results sug- 
gest that the investigation of the SC phase under the 
pressure would be interesting, in order to explore the SC 
phase in the regime where AFM order does no longer 
exist. 
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ful discussions. This study has been supported by the 
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and the NCCR MaNEP. Numerical computation was car- 
ried out at the Yukawa Institute Computer Facility. 



It has been reported that the Knight shift remains con- 
stant below T c for any field direction. 21,22 This result 
looks puzzling in view of calculations which suggest the 
decrease of spin susceptibility for fields in the a6-plane to 
the half of its normal state value. 8 ~ 10, 13 Here, we point 
out that this discrepancy can be resolved by taking into 
account the AFM order. The uniform spin susceptibil- 
ity of the normal and SC states at T = is shown in 
Fig. 4, assuming T c -C a. No correlation effects have 
been taken into account here. For fields H _L Hq the nor- 
mal state and SC state susceptibility merge for increas- 
ing staggered moment. On the other hand, for H \\ Hq 
the behavior is opposite. Assuming that the anisotropy 
energy is sufficiently small, the condition H _L Iiq is gen- 
erally favored. Together with other effects such as vortex 
scattering and the formation of a helical SC phase 12 the 
influence of AFM order could eventually account for the 
experimental results. 21 ' 22 Note that Fig. 4 is obtained 
without taking into account the canting of AFM moment 
in contrast to Ref. 15 where the AFM coupling to the 
local moment is assumed to explain a similar property. 
Another mechanism to enhance the spin susceptibility 
below T c has been proposed by Fujimoto, which is based 
on a strong particle- hole asymmetry. 14 However the j3- 
band in our model eq. (1) does not satisfy the necessary 
conditions. If the AFM order is the main cause for the 
behavior of the spin susceptibility for in-plane fields, a 
distinct change should be observed when the AFM order 
is suppressed by pressure. 

In summary, we have examined various aspects of the 
pairing state in CePtaSi based on a Hubbard model in- 
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